Magnetic Excitations of Stripes Near a Quantum Critical Point 
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We calculate the dynamical spin structure factor of spin waves for weakly coupled stripes. At 
low energy, the spin wave cone intensity is strongly peaked on the inner branches. As energy is 
increased, there is a saddlepoint followed by a square-shaped continuum rotated 45° from the low 
energy peaks. This is reminiscent of recent high energy neutron scattering data on the cuprates. 
The similarity at high energy between this semiclassical treatment and quantum fluctuations in 
spin ladders may be attributed to the proximity of a quantum critical point with a small critical 
exponent n. 
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Strongly correlated electronic systems often exhibit 
some evidence of local electronic inhomogeneity, appear- 
ing in such diverse probes as neutron scattering, STM, 
/isR, and NMR, among others. 0, 0, Q In the high tem- 
perature superconductors, stripe structures are one pos- 
sible microscopic realization of local electronic inhomo- 
geneity. Stripes are electronic states which spontaneously 
break the rotational symmetry of the host crystal, the 
most ordered example being an interleaved, unidirec- 
tional modulation of both spin and charge density, as 
shown in Fig. ^ Recent neutron scattering experiments 
exploring the spin excitations in cuprates point to a uni- 
versal high energy magnetic response in both YBCO and 
(stripe-ordered) LB C 00 HQ, and it has been suggested 
that the universal high energy behavior may be connected 
to the superconductivity in these materials. Q At low 
frequency, neutron scattering reveals four incommensu- 
rate spin peaks which disperse inward toward the (tt, tt) 
point, converging at a resonance peak at intermediate 
energy. This is followed by a high energy square shaped 
continuum in which the corners are rotated 45° away 
from the direction of the low energy peaks. [113] The high 
energy response in particular has been attributed to the 
quantum excitations of spin ladders. 0H, 0,0 

We find that the experimental results in LBCO at all 
energies are also consistent with semiclassical spin wave 
excitations of weakly coupled stripes, i.e. for weak spin 
coupling across the charged domain walls in the spin pat- 
tern. We further suggest that the reason the quantum ex- 
citations of 2-leg ladders and the semiclassical spin waves 
studied here have such similar behavior at high energy is 
due to the proximity of a quantum critical point (QCP) 
with small critical exponent r\. On both the ordered and 
disordered side of a QCP, above a certain characteristic 
frequency the response is quantum critical, but for small 
77, the quantum critical behavior can look very much like 
the Goldstone behavior on the ordered side.||| Calcula- 
tions of quantum spin fluctuations due to high energy 
ladder behavior jfj,|fj,l3| are exploring this quantum criti- 
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FIG. 1: (a) Site-centered stripes of spacing 4 viewed as weakly 
coupled 3-leg ladders, (b) Bond-centered stripes of spacing 4 
viewed as weakly coupled 4-leg ladders. The coupling within 
each ladder is J a , and the coupling between ladders is Jj,. 



cal regime at high energy, and can resemble much of the 
semiclassical behavior reported here. 

In this Letter, we consider fully ordered static spin 
stripes, which are arrays of antiphase domain walls in an 
otherwise antiferromagnetic texture. We are interested 
solely in the response of the spin degrees of freedom and 
neglect the dynamics of the charge density which must 
peak on every domain wall. These charge degrees of free- 
dom affect the spin degrees of freedom through a spatially 
modulated effective exchange integral. The undoped sys- 
tem was shown to be quantitatively well described b 
the semiclassical spin waves of a 2D antiferromagnet. 
When stripes form, the doping is topological, that is, 
doped holes create line dislocations in the antiferromag- 
netism, but in between the defects the antiferromagnet is 
more or less intact. Within this framework, it is reason- 
able to expect a semiclassical treatment to be applicable 
in a range of doping, as long as the ground state remains 
ordered. We take the coupling within an antiferromag- 
netic patch to be close to its full undoped value, and 
consider the effective spin coupling across a charged do- 
main wall to be reduced, as illustrated by the dotted lines 
in Fig. ^ This corresponds to a Heisenberg model with 
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modulated exchange integral on a square lattice, where 
each site represents a copper atom in the copper-oxygen 
plane: 



H = \ E J "'( 



S r • S r ') 



(1) 



<r,r') 



where J rr < is the exchange coupling. Nearest neighbor 
couplings are positive J r , r > — J a > within each anti- 
ferromagnetic patch. Couplings across domain walls de- 
pend on whether stripes are site- or bond-centered, as 
explained below. We work in units where ft, = 1. 

"Site-centered" stripes have domain walls which are 
centered on the sites of the square lattice (i.e. on the 
copper sites), leading to an antiferromagnetic effective 
coupling Jr !r / = Jb > across the domain walls. For 
weak coupling across the domain walls Jb <SC J a , the sys- 
tem is close t o the regime of coupled 3-leg ladders, as 
shown in Fig. 1(a) "Bond-centered" stripes have do- 
main walls which are centered between the sites, lead- 
ing to a ferromagnetic effective coupling J r r > = Jb < 
across the domain walls. 01 F° r weak coupling across the 
domain walls \Jb\ <C J a , the system is close to coupled 
4- leg ladders. jllj as shown in Fig. 1(b) When Jb is suf- 
ficiently large, both cases presumably have a long-range 
ordered ground state^ with low energy Goldstone be- 
havior, but there is a quantum critical point at small Jb 
beyond which the physics scales to that of decoupled lad- 
ders. While the ground state of coupled odd-leg s = 1/2 
ladders is ordered with any finite coupling and therefore 
J"' 1 * = 0, coupled even-leg ladders have a finite value for 
the quantum critical point, |J£ rit | > 0. 

We use semiclassical linearized spin wave theory and 
Holstein-Primakoff bosons, a standard procedure de- 
scribed elsewhere p^l Hil llij . in order to calculate the 
spin wave excitation spectrum and the zero-temperature 
dynamical structure factor, 

s(k,cu) = j2 E K/|fiW>i 2 *(w-w/) (2) 

/ i=x,y,z 

which is related to the expected neutron scattering in- 
tensity. 

Fig. shows the intensity of the dynamical structure 
factor for the site-centered stripes of Fig. 1(a) as a func- 
tion of frequency for weak coupling strength Jb = 0.05 J a . 
Fig. |31 shows the results for the bond-centered stripes of 
Fig. 1(b) for weak coupling strength J& = — 0.09J a . Both 
coupling strengths have ordered ground states, as dis- 
cussed below. The results are rather similar for site- and 
bond-centered stripes, Q with the main difference being 
that while the satellite peaks are too weak to be visible 
at low energy in the site-centered case, they are visible 
(although faint) at low energy in the bond-centered case. 
At low energies, four strong incommensurate peaks are 
visible, which disperse inward toward {%, it) as energy 
is increased. Although a spin wave cone must emanate 
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FIG. 2: (Color online) Constant energy cuts, integrated 
within an energy window of ±0.05J a S for vertical, site- 
centered stripes of spacing 4 at Jb = 0.05J a in the magnetic 
Brillouin zone. The energy E is in units of J a S. Results are 
shown for twinned stripes. 



from each magnetic reciprocal lattice vector due to Gold- 
stone's theorem, the intensity is not necessarily uniform. 
For weak spin coupling across the charged domain wall, 
\Jb\ *C | J 1, we find that the intensity is strongly peaked 
on the inner branches emanating from (jr, it ± 7r/4), and 
the twinned IC peaks at (tt ± 77/4 ,71") , that is, the side 
of the cones that is nearest (tt. tt) . |lq] This situation is 
reversed for strong \Jb\ ^> \J a \, where the intensity in 
the spin wave cones is strongest on the outer branches 
emanating from the same points. One mystery about 
the low energy neutron scattering results in cuprates has 
been that spin wave cones are not observed, but rather, 
the intensity disperses toward (7r,7r). We have shown 
here that this is consistent with semiclassical spin waves 
of coupled stripes. 

Upon increasing the energy, there is a resonance peak 
apparent at E les = 0.63 J a S and (tt, it) in Fig. El which is 
a saddlepoint at which the integrated intensity S(ui) has 
a maximum (see Fig. 0J. In Fig. |3 the resonance is at 
E les = 0.75J a S. The energy scale of this resonance in- 
creases with increasing \Jb/J a \- Using E rcs — 50 — 60meV 
from Ref. Q, this yields J a = 160 - 190meV in Fig. 
and J a = 130 — 160meV in Fig. |SJ This is reasonable 
given a value of J a — 140meV in undoped La2CuO4 0], 
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FIG. 3: (Color online) Constant energy cuts, integrated 
within an energy window of ±0.05 J a S for vertical, bond- 
centered stripes of spacing 4 at Jt, = — 0.09J a in the magnetic 
Brillouin zone. The energy E is in units of J a S. Results are 
shown for twinned stripes. 




FIG. 4: S(ui) for site-centered stripes of spacing 4 at Jt — 
0.05J a over the magnetic Brillouin zone. The frequency lu is 
in units of J a S. The broadening Acu is of ±0.05J a S. 



and it is consistent with other theoretical estimates on 
LBCO.jfl Q For weak enough coupling J^, the disper- 
sion above the saddlepoint is highly anisotropic, giving 
rise to a rotated square-shaped continuum above the res- 
onance peak, much like what is seen in LBCO.jlJ. 

Fig. 2| shows the single magnon contribution to the 
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FIG. 5: Crossovers for coupled ladders as a function of inter- 
ladder coupling Jt for both odd-leg ladders (a) and even-leg 
ladders (b). 



momentum-integrated structure factor S(u>) for site- 
centered stripes of Fig. O and reveals the effect of the 
saddlepoint in the acoustic band. The broad peak in 
S(uj) at E les = 0.63 J a S is due to the saddlepoint. Quali- 
tatively, the results at high energy are similar to the one- 
triplon calculations based on coupled 2-leg ladders. [HQ 
Contributions from the first optical band can be seen 
above E = 2.2J a S. 

We find that we can describe both the low energy 
and the high energy neutron scattering data on LBCO0] 
within the semiclassical framework of spin waves, for 
weak spin coupling across the charge domain walls. This 
brings up two questions: 1. Why are the semiclassical re- 
sults so similar to the quantum critical behavior in other 
models of stripes, 0] and 2. Why is the high energy 
response of YBCO so similar to that of stripe-ordered 
LBCO, when YBCO does not show evidence of long range 
magnetic order? Concerning the first question, the sim- 
ilarity may be due to the proximity to a quantum crit- 
ical point with small critical exponent rj. || The ele- 
mentary excitations of the system fundamentally change 
in character across the quantum critical point, and for 
rj > the structure factors in the ordered and quantum 
critical regimes have analytically distinct forms. How- 
ever, the spectrum is distributed in a similar way for 
77 > 0. In fact, as i] — > 0, the two cases become indis- 
tinguishable. We are considering a 2D quantum transi- 
tion to distinct 3-leg or 4-leg ladders, and the universal- 
ity class is the classical (2 + l)-dimensional Heisenberg 
model. Since then the critical exponent 77 = 0.037 is 
small, for all practical purposes there may be little 
distinction between the classical and quantum critical 
cases. However, sufficient resolution in lineshapes can 
in principle distinguish: whereas weakly interacting spin 
waves produce a Lorentzian lines hape, q uantum critical- 
ity produces a power law cusp, [a, llTJ, l2f| Concerning the 
second question, the gap in YBCO indicates that it must 
be on the disordered side of the QCP. However, the fact 
that YBCO's high energy response is so similar to the 
Goldstone modes calculated here is likely an indication 
that YBCO is close to the QCP. 



Fig. 5(a) shows the crossover energy scales for coupled 
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odd- leg ladders as a function of interladder coupling 
In this case, the quantum critical point is at J" 1 ' = 0, 
where the system breaks up into independent odd-leg 
ladders for S = 1/2. At = 0, the low frequency be- 
havior is controlled by the c — 1 Wess-Zumino-Wittcn 
(WZW) model, and the high frequency behavior scales 
to that of weakly interacting Goldstone modes. 01 Away 
from the QCP, the ground state is ordered and the low- 
est energy response is due to Goldstone modes. Quan- 
tum critical behavior obtains above a lower crossover 
(Jfc)^ 2 ,!!^! and the high energy Goldstone be- 
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havior of decoupled odd-leg ladders is recovered above 
a higher crossover energy scale luq oc j oc 



-ir, 



where the effective width of ladders W e f f increases as Jf, 
increases. The two energy scales Wq C and Wq approach 
each other with increasing Jb/J a , so that beyond a cer- 
tain point the WZW behavior is squeezed out entirely, [l^ 
and in this region our approach is well-justified. More- 
over, because the critical exponent r\ = 0.037 is small, 
even the quantum critical WZW region can resemble the 
Goldstone behavior. 81 



Fig. 5(b) shows the gap A and the quantum critical 
energy scale Wq C for coupled even-leg ladders as a func- 
tion of J&.[2l| In this case, the quantum phase transi- 
tion is at finite coupling J"' 1 * > 0. To find the quan- 
tum critical point for coupled 4-leg ladders, we use the 
stochastic series expansion quantum Monte Carlo (QMC) 
method. [13] Using finite size scaling on the spin stiffness 
and the Binder ratio, we find that J£ nt = 0.076(3) J a for 
weakly coupled 4-leg spin- 1/2 ladders, consistent with 
Ref. [23]. (For weakly coupled 2-leg ladders, J" 1 * « 
0.25J a — O.3J O .0,|3) Ovrr spin-wave calculations pre- 
sented in Fig. |3] correspond to weakly-coupled 4-leg lad- 
ders, close to but on the ordered side of the QCP. Using 
QMC on this configuration for systems up to size 64 x 32, 
we find that the sublattice magnetization is m z = 0.143 
(47% of 0.307 for the 2D antiferromagnet) , confirming 
that the system is ordered at this coupling. 

Quenched disorder in the form of dopant atoms intro- 
duces disorder in the ladder widths, positions, and cou- 
plings. Weak disorder is irrelevant at the critical point 
of a (2 + l)-dimensional Heisenberg model, and although 
disorder can lower the transition temperature (a finite 
ordering temperature is possible with any weak c-axis 
coupling 9]), the critical exponent 77 is unchanged. Dop- 
ing may even induce a QCP through the introduction of 
disorder, doping the ladders, or by inducing dynamical 
stripes with fluctuating charge. 

In conclusion, we have studied the magnetic excita- 
tions of ordered stripes in the proximity of a quantum 
critical point. For weak spin coupling across the charged 
domain walls, the low energy spin wave cones have weight 
which disperses inward toward (tt, it), in agreement with 
experiment. At higher energies, the semiclassical excita- 
tions develop saddlepoints with steep dispersions. The 
saddlepoint acts like a "resonance peak" in that there 



is extra intensity due to the saddlepoint, and also due 
to stripe twinning. (However, our model does not ad- 
dress the observed increase in intensity of the resonance 
peak as superconductivity onsets.) At energies above the 
saddlepoint, there is a slowly dispersing square-shaped 
continuum. The proximity to a quantum critical point 
(either to decoupled 3-leg ladders or to weakly-coupled 
4-leg ladders) with small critical exponent rj implies that 
quantum critical behavior (and therefore quantum lad- 
der calculations at high energy) can strongly resemble 
the semiclassical calculations presented here. 
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